Relationships among coefficients in deterministic and stochastic transient diffusion 

Z. Kaufmann 

Institute for Solid State Physics, Edtvos University, Muzeum krt. 6-8, H-1088 Budapest, Hungary 



ON 
On 



Oh 
00 



> 
O 

o 

On 
O 

t> 

ON 
I 

O 

a 

43 

o 



Systems are studied in which transport is possible due to large extension with open boundaries in 
certain directions but the particles responsible for transport can disappear from it by leaving it in 
other directions, by chemical reaction or by adsorption. The connection of the total escape rate, the 
rate of the disappearance and the diffusion constant is investigated. It leads to the observation that 
the diffusion coefficient defined by {x 2 ) is in general different from the one present in the effective 
Fokker-Planck equation. The result makes it possible to generalize the Gaspard-Nicolis formula for 
deterministic systems to this transient case. 



Randomness and diffusion are common features of ex- 
tended stochastic and chaotic systems IQg^^ggQgg^ . 
Among deterministic systems diffusion has been much 
studied in the Lorentz gas ||[| . As more simple models 
proper ID maps |^,^|,0J§| and a 2D map || has been in- 
troduced, which were built as chains of maps. Ref. fic| l 
has shown a relationship between the diffusion coefficient 
and microscopic quantities, namely the Liapunov expo- 
nent and the Kolmogorov-Sinai entropy referring to the 
repeller. This was later generalized to other transport 
coefficients |fl| and case of small external field . The 
Liapunov exponent was independently calculated for the 
random Lorentz gas |l3| . 

There are systems in which particles can escape in di- 
rections transversal to the extension of the system raising 
interestic problems |H] in the field of transient chaos |l5| ] 
in particular in the critical case fl(^Jl7|] . A simple example 
is a channel in a mesoscopic system modeled by a strip 
of Lorentz gas with open side boundaries. Another ex- 
ample is the Troll-Smilansky model for chaotic scattering 
consisting of a ID infinite periodic array of soft potential 
. !-h ! valleys, which is a system related to a model of ionization 
|p^ , ^9| . Further example is the motion in an infinite set 
of resonances in the phase space of Hamiltonian systems 
when one considers the particle to escape when it leaves 
the set of chosen resonances |2Q] . Diffusion has also been 
investigated on a chaotic saddle in a model for the in- 
teraction of a particle with electrostatic wavepacket [^l] . 
Particles can also be lost from the point of view of dif- 
fusion by absorption, chemical reaction p2| or by other 
ways. Particles in such transiently chaotic or stochastic 
systems can diffuse in the extended direction or direc- 
tions for some time and then escape either through the 
ends in the extended direction (if the system is finite) 
or in the other directions or other ways. Therefore the 
average period for which they take part in transport is 
finite, and remains finite even in the limit when the size 
of the system in the extended direction goes to infinity. 

The aim of the present paper is to generalize the 
Gaspard-Nicolis formula |nj to the case of the above de- 
scribed transient diffusion in deterministic systems. For 



this purpose it is studied how the total escape rate sepa- 
rates to terms related to the extended and the transver- 
sal direction. This investigation is made in a general way 
leading to interesting results for both deterministic and 
stochastic systems. For simplicity the system is assumed 
to possess one extended direction with a discrete transla- 
tional symmetry, and in most of the considerations here, 
an inversion symmetry that reverses the extended direc- 
tion. For sake of convenience the primitive cells with re- 
spect to the translational symmetry shall be labeled by 
a discrete variable x that is monotonous in the extended 
direction. The rest of variables specifying the state of the 
particle shall be assembled in y. It is convenient to show 
the choice of x and y in case of the strip of Lorentz gas. 
Here, as in general, it is easier to study the system in 
discrete instead of continuous time. Taking the surfaces 
of the disks as Poincare surface the state of the particle 
on it can be given by the coordinates x, q, a and (3. By x 
we mean the ordered label of the periods of the structure. 
q is a label of the disk inside one period, a is the angle 
of position on the disk and (3 is the angle of reflection. 
Then y corresponds to (g,a,/3). 

The general evolution equation for the probability den- 
sity </>t(x, y) of the particle can be written as 



fa+i(x,y) = 



J y> 



u 3,v,v 



(f>t(x-j,y') , 



(1) 



where, if y contains continuous variables one can consider 
the sums as integrals over the continuous components 
or one can use coarse graining with arbitrary precision. 
Here (fit = if the argument falls outside the region of 
the system. The maximal jump J can be assumed to be 
finite, or the transition probability w^ v ^i to decay fast 
in j. The translational symmetry is implied in the form 
of while an inversion symmetry can be written as 
w j,y,y' = w -j,T y ,Ty', where T 2 y = y for every y. 

Two representative classes of such systems can one 
keep in mind here. The first is a 2D random walk in 
a strip, for simplicity assuming no memory, for which 
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0,y) = 



K 

E 

— J k=-K 



W jk cj> t {x-j,y-k) (2) 



applies. This can be considered as a rough descrip- 
tion of the Lorentz gas strip. The inversion symme- 
try can be a point inversion or a line inversion symme- 
try. To give an account of correlation between tran- 
sitions ID walks with memory are taken as a second 
group of examples. Then ([!]) can be used conceiving 
y as the memory containing, say, n number of past 
steps ji,j 2 , .-.j n and% +l! „- = P(j n+1 \j n , j„_i, . . . ji) 
is the conditional probability of the next step. The 
inversion symmetry implies P(j n +i\jn, jn-i, ■ ■ ■ h) = 

P (-jn+l\ ~ j„, -j„_l, . . . - Jl). 

For general considerations we shall return to (Q) taken 
without restriction to ID. For the diffusion process the 
long time behavior of the system is important. That is 
governed by the leading eigenfunction of the right hand 
side of (|l|), i. e. the solution 4> t for which cf> t+1 = e^ K (f> t 
and the total escape rate k is minimal. This 4>t shall 
be called the asymptotic state. The boundary of the 
system should also show the symmetry, so the region of 
the system should be defined by independent conditions 
in x and y (x £ R x and y £ Ry)- If the particle from a 
point (xo,yo) jumps to a point for which y Ry 

the particle shall be considered to escape in y direction, 
while in case y € R y but x ^ R x it escapes in x direction. 

To find the relation of the escape rates and the diffusion 
coefficient is simple in certain cases but is problematic in 
general. It is enlightening to study the simple cases first. 
These are the ones in which the asymptotic solution sep- 
arates as a product <p(x,y) — tp(x)co(y). This happens if 
the transition probability matrix is a sum of diadic prod- 



ucts Wjyy, =J2s u 



(«)„,« 



? v yy anc ^ ^ e P ar ^ a ^ diffusion equa- 



tions in x direction have a common leading eigenfunction, 

i. e. J2j u^ip{ x ~ j) = fJ-aip(%) such that for every s /x s is 
the maximal eigenvalue. This property shall be denoted 
by S x . It is easy to see that in this case the eigenvalue 
equation in y direction v yy l fj, s u(y') = e~ K uj(y) de- 
termines uj and e~ K . An example for this case among the 
2D walks (H) is a walk on a square lattice [1,L] £5 [1, M] 
with 

{WjkVA -i=\ r ) = ( r |o(010) + (0 )o(101) (3) 
\p p \0 \p 



whose solution reads 
2£ 



Another possibility is when the eigenmode is common 
in y direction, i. e. Y, y > v y y '^{v') = v s^{y) (property S y ). 
Then JZj a u^v s ip(x — j) = e~ K ip(x) determines ip and 
e~ K . When the properties S x and S v are both satisfied 
then e~ K = J2 S ^sV s - 

To obtain an evolution equation of the distribution on 
large scale in x direction in either of the above cases with 



4>(x,y) = ifj(x)ui(y) one has to study the evolution of %[>. 
It is convenient to choose the normalization uj(y) = 1, 
since then ip(x) — J2 y ^( x ^ V)- Using this equation and 
(Q) one obtains 

i>t+i{x) =^2<w j il) t (x-j) , Wj =^2w jt y iyl u(y') . (4) 
j yy' 

Obviously the probability that a particle at xq = x — j 
does not escape in y direction in the next step is e~ Ky — 
Sj^j = J2jyy' w jyy ,LU (y')> which is independent of x . 
One can separate this escape with substitution ipt{x) = 
e^^gtix), which yields g t +i(x) = J2j ^j^" 9t(x - j). 
Since this equation describes a random walk and the sys- 
tem is extended in x direction an effective Fokker-Planck 
equation g t +i(x) — gt{x) + Dppg'tix) is valid on large 
scales. Returning to ip it takes the form 



^t+i{x) = ( 1 + Dfp-^ ) e- K "Vt 0) 



(5) 



Its dominant solution in case of a channel of length L is 



e~ Kt sin 



in (4) 



(6) 



and k x = (n 2 /L 2 )D F p + 0(L~ 3 ), if defined such, 
that e _Kx is the conditional probability that a point does 
not escape in x direction if it does not escape in y direc- 
tion. Note, that k v may also depend on L. 

The general case, when cj)(x,y) ^ ip(x)Lu(y) is more 
complicated. The main point is, that n v becomes depen- 
dent on x and L, the length of the system in x direction. 
However, in case when L is much larger than the size 
in other directions and there is an inversion symmetry 
it shall be shown, that apart from a vicinity of the ends 
of the channel the deviation of k v from a value Ky°°^ 
is proportional to f"(x)/ f(x). Here f{x) is introduced 
analogously to ip(x) as f(x) = J2 y 4>{ x , y) an d Ky 00 ^ is the 
value of k v for the homogeneous solution in case L — oo. 
This makes it possible to write down a proper effective 
Fokker-Planck equation. It is suitable to choose a seg- 
ment of the system separated by x — const planes, such 
that its size is still much larger than the transversal size, 
but much smaller than L. Such a segment feels values 
f(x\), f{x2) of / at its ends with a common exponen- 
tial decay e~ Kt . It is assumed that the diffusion in the 
system mixes the contribution of sites with different y co- 
ordinates. So systems are excluded in which the particles 
from one site can not fill the whole system, thereby differ- 
ent initial distributions can lead to different asymptotic 
states. With this assumption the effect of the transver- 
sal distribution <p{x, y)/f(x) decays fast from the ends 
of the segment towards its inside. On the other side, in 
the asymptotic state the transversal distribution at the 
ends of the segment is almost identical to the one in the 
middle. Therefore one can expect the distribution in the 
inside of the segment is determined by the values f(xi), 
f(x2) and the decay rate k. 
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Thereby the distribution can have three free param- 
eters. Since the evolution equation (Q) is linear in the 
distribution, the asymptotic distributions <fi and / can 
be assumed to have linear combination form 

<t>(x,y) = [(/) {oc) (x - x ,y) + ci(f> (1) (x - x ,y) 

+c 2 ^(x-x ,y)}e- Kt , (7) 
f(x) = [1 + Cl /«(x - xq) + c 2 f {2 \x - x )]e- Kt , (8) 

where fw(x) = < - fc - > (x, y) for k — 1,2 and xq is the 



center of the segment. Here </>(°°) (x, y) 



(oo) 



(y) is the 



oo that is independent of x. 



solution in the limit L 

The corresponding escape rate shall be denoted by Ky 
cj)^ is the asymptotic solution of ([!]) with antisymmet- 
ric boundary conditions f{x\) — —f(xi). This term is 
responsible for current through the middle of the seg- 
ment. Starting with a </>(°°) alone and symmetric bound- 
ary conditions f{x\) = /(a^) = ae~ Kt one observes in 
general that <j> in the middle decays faster or slower than 
on the boundary depending on the sign of k— k^°°\ In the 
asymptotic state this leads to a hump- or vale-shape term 
f( 2 ^ in /, which corresponds to some 4>( 2 \ The partial 
distributions can be approximated as (x—xq) = x—xq 
and f( 2 '(x— Xq) — fx— xo) 2 — b in a vicinity of the middle 
of the segment if 4> , k = 1,2 are properly normalized. 
This means, that the general solution is 

4>(x, y) = f(x )(/)^ (x-Xq, y) 
+ f'(x )4>W(x-x ,y) 

+ ^^{4> (2) {x - x , y) + b^\x - x , y)) (9) 

in a vicinity of the point Xo- The local rate of escape 
Ky(xo) in y direction in the middle can be calculated as 



1 ,-«,(«) - E ^ 
f(xo) ' 



(10) 

E ( x o) = X! 0( X O' v)-^2 w jyy ><f>(x 0) y') , (11) 
v yy' 

where E(xq) is the flow of escape in y direction at Xo- 
Then clearly 

E(x ) = E oo f(x ) + E 1 f'(x )+ {E2+ ^ Eoo) f"(x ), (12) 

where E&, k = oo, 1, 2 are constants characteristic of (j>^ , 
respectively. It can be shown that E\ = 0. To see this 
one can separate Eq. (|l|) to terms related to transition 
from inside of the chosen segment and from outside. Con- 
ceiving (f>t(x,y) as a vector $ = where any value 
of the index i corresponds to a point in (x,y) space in- 
side the segment and assuming the asymptotic state (i. e. 
4>t+i = e~ K t ) one obtains e _K $. i = J2j Wij'&j+Xi- Here 
w describes the transitions inside, x the transitions from 
outside, and w, \ can be constructed using w and values 
of <j> outside but near the boundary of the segment. The 



solution for $ is <f> = (e~ K I — w) x i where / is the 
unit matrix. The boundary conditions are antisymmet- 
ric for the symmetry transformation T in the state ft 1 ' , 
thereby x is a l so antisymmetric, w and w are symmetric 
for T. Consequently $W and the corresponding <jP^ are 
antisymmetric. Therefore ( |TT|) yields E\=Q. Using jl0| ) 
and (0) one obtains 



-(oo) 



fix) 



+ G(L- 3 ) 



(13) 



with a suitable constant ry, since for large L one expects 
/" = 0(L~ 2 ). This is not valid near the ends of the chan- 
nel where the middle of the segment can not be placed. 
Then the analog of the effective Fokker-Planck equation 
(ft) becomes of the from 



ft(x) , 



or, with substitution of ( |l5| ) 



/■/+i<,'\> = < " [l + (D F p-v)^)ft(x) 



(14) 



(15) 



It is plausible to assume that a distribution corresponding 
to (|9J) sets in earlier in time than the asymptotic state 
in the extended direction. Then (|lj,|l^) are also valid 
for general /(x) distributions still not in the asymptotic 
state. ( |l5| ) shall be applied later for such a case, but here 
the asymptotic solution is important. It is given by 



ft(x) 







Hi' 









(16) 



with k — Ky — log[l — {Dpp — r])a / L ■], which is valid 
inside the channel with a deviation at its ends, such that 
ft(x) reaches zero at the ends (x = and x = L), while 
( Jig ) takes zero value at £ = (1 — n/a)L/2. In the limit 
L — > oo the value of £ becomes constant, since the neigh- 
borhood of the channel behaves in the same way in this 
limit. Therefore a = n + C(L _1 ). Using @ and @ 
one can notice, that k v (x) is constant apart from the 
vicinity of the ends of the channel. Its value is 



XL) 



7T 
L2 



v + o(l- 3 ) , 



(17) 



where the notation emphasizes its dependence on L. Us- 
ing this equation and the ones for k and a one obtains the 
relation of the escape rates with the diffusion coefficient 



tD 



FP 



0(L~ 



(18) 



Note, that (tt 2 /L 2 )Dfp does not give k x in general, con- 
trary to (§). Following |l(J in case of deterministic pro- 
cess the total escape rate can be related to the Liapunov 
exponents and the Kolmogorov- Sinai entropy of the re- 
peller, namely, k = J2\ i>0 ^i — hxs f^^fj- This yields 
the generalization of the Gaspard-Nicolis formula PQ] 



3 



A,>0 



h KS = J1 d fp + 4 L) + 0(L- 3 ) . (19) 



One more important question is whether Dpp de- 
fined as the coefficient in (|l4|) is equal to the one de- 
fined by the mean square deviation of a: as J ft (x) (x — 
xo) 2 dxj J ft{x)dx oc 2D ms dt, starting from a state con- 
centrated in a vicinity of xq. Introducing gt{x) — 

e K " ]t ft{x) in Eq. ( |l5| ) one can eliminate the factor con- 
taining Then the equation becomes an evolution 
equation of a diffusion process whose diffusion coefficient 



D 



msd 



= Dpp - 7] 



(20) 



is clearly equal to the diffusion coefficient defined for 
ft(x) by mean square deviation. So we can see Dpp 
and D ms d are in general different. Comparing (|l9|) to 



the Gaspard-Nicolis formula |10[ one can notice that k 



(L) 



has appeared as an additional term, but using (17) and 
( p0| ) one can get another form of (|l^) in which D ms d and 

are present instead of Dpp and k 



(L) 



E 

Ai>0 



A, - h KS = j^D msd + 4°°) + 0(L- 3 ) 



(21) 



Numerical calculations have been made to test the va- 
lidity of the statements for k v (x) and Eqs. (|l4|,|o|) in two 
concrete models. 

A, The first model is a 2D random walk on a square 
lattice [1,L] <g> [1,M] with 



{Wjk}-{,-i = I 9 r q 




(22) 



and p= 0.1, 9 = 0.2, r = 0.4. 

B, The second model is a ID random walk with two 
step memory such that P(j t +i \jt, jt-i) = Rjdt-iQh+ijt 
and R +1 = 0.45, = 0.9, Q +1 = 4/9, Q_ x = 5/9. To 
ensure symmetry Q and R depend only on the product 
in their subscript. Q + i + Q_i = 1, therefore R describes 
the probability that the particle does not escape in y 
direction and Q describes the relative probability of the 
steps jt+i- A possibility was given to rest for one step 
with a probability g = 0.01 without changing jt,jt-i- 

In model A K y (x) was found to decay exponentially 

to Ky coming from each of the ends towards the in- 
side of the region. In model B K v (x) is constant for 
x = 2, 3, ... L — 1 and has different values only at the 
endpoints. In both models k and k v {x) have been mea- 
sured. Dpp has been calculated from them by ( |Is| ) 

for large L and r\ from the L-dependence of n y L ^ by 
(|l7j). D ms d has been determined independently measur- 
ing ((x — xo) 2 } for a well concentrated initial distribution. 
The results satisfy (|(j|) up to the expected precision (5 
digits). In case of model A Dpp has also been calculated 
by i(x) = D F p(f(x) - f(x + l))e~ K y( x \ where i(x) is the 



current between sites at x and x + 1 . The results of this 
match the ones by ([l8]), thereby confirming the validity 
of (HJ) and (|H) , which was used in (|l9 ) . 
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